Long-time signatures of short-time dynamics in decaying quantum-chaotic systems 
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We analyze the decay of classically chaotic quantum systems in the presence of fast ballistic escape 
routes on the Ehrenfest time scale. For a continuous excitation process, the form factor of the decay 
cross section deviates from the universal random-matrix result on the Heisenberg time scale, i.e. for 
times much larger than the time for ballistic escape. We derive an exact analytical description and 
compare our results with numerical simulations for a dynamical model. 



Properties of complex quantum systems are mostly 
characterized by a high degree of universality which 
is rooted in the random interference of many partial 
waves 0. A prototypical example arc the cross section 
fluctuations which arise from the random interference be- 
tween the decay modes of an excited quantum system. 
In the case of a large number of decay channels, these 
fluctuations are free of all characteristics except for the 
classical decay rate (Ericson fluctuations) Mi Q . The 
main theoretical tool to describe such universal proper- 
ties is random-matrix theory (RMT) [^|^. Semiclassical 
theories reveal limits to this description which are rooted 
in direct processes and hence best characterized in the 
time domain. Most notably, if the complex quantum dy- 
namics is due to a chaotic classical limit, the random 
wave interference is only established after a finite time, 
the so-called Ehrenfest time 0. For shorter times, the 
quantum dynamics follows the classical dynamics quasi- 
detcrministically. Over the past years, a large number of 
studies have explored the consequences of a finite Ehren- 
fest time H, 0, 0| I but almost exclusively these stud- 
ies were concerned with stationary properties. 

The purpose of this paper is to develop a dynamical 
theory for a physical observable - the form factor of the 
decay cross section of an excited quantum system - in 
which non-universal corrections can be detected directly 
in the time domain. The dynamical theory is based on 
a stroboscopic description which extends the previous 
analysis of stationar y p roperties in systems with a fi- 
nite Ehrenfest time [iSj- In the absence of direct 
decay processes, our model exhibits the typical univer- 
sal fluctuations found in autonomous models of quantum 
decay and described by RMT 0, . We incorporate 
quasi-deterministic direct decay by applying a mapping 
formalism which has earlier been used to describe Fano 
resonances due to direct decay processes 0, 0, . 
It turns out that the deviations in the form factor from 
the universal behavior persist up to times comparable 
to the Heisenberg time, hence far beyond the regime of 
quasi-deterministic decay. 

Molecular photo-dissociation and atomic auto- 
ionization are examples of half collision processes (see 



Ref. |20j and references therein), i.e. the coherent 
dynamics of quantum systems which have been excited 
by external means {e.g. the absorption of a photon) 
and decay subsequently via energetically accessible open 
channels. Within the half-collision description, each 
decay amplitude is given by the overlap {^a\oi) of the 
square integrablc initial state \a) and the scattering 
state |4'a) associated to channel a. The form factor 
C{t) is defined as the Fourier transform of the auto- 
correlation function of the total decay cross section 
a{E) = K^al*^)!^- It can also be obtained from the 
Fourier transform of the cross section itself. 



C{t)^^\ait)\\ d{t) = 



L/2 

dE e 

-L/2 



-27ri Et 



a{E). (1) 



Here, we measure time in units of the Heisenberg time 
tn = h/ IS., and energy in units of the mean level spacing 
A {h is Planck's constant). Typically, only a finite energy 
range, classically small but quantum mechanically large 
[L 1), enters this formula. 

In previous works Q, 0, 0, 0, 0|, the random- 
matrix description of half-collision processes has been 
set up for autonomous systems which are described by 
an effective Hamiltonian jj. In the case of time rever- 
sal invariance, the autocorrelation function then follows 
from the Verbaarschot-Weidenmiiller-Zirnbauer (VWZ) 
integral and for systems without direct decay the 
form factor is given by 
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Stroboscopic model. We base our dynamical theory 
of half-collision processes on a time-periodic quantum 
map description which has been previously used to de- 
scribe transport in systems with a finite Ehrenfest time, 
^Ehr llJ, UlA I2JI ■ Because of time-reversal invari- 
ance of the decay dynamics, the Floquet operator F is 
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symmetric, F = F^, which allows to decompose F into 

F = Fout Fin, with [Finf = Font 0. To introduce decay, 

we specify an TV-dimensional subspace (spanned by a col- 
umn vectors of the matrix p) within the M-dimensional 
Hilbert space of F. That subspace provides the interface 
between the closed system and the decay channels. For 
baUistic decay (ideal coupling) we define Q = 1—pp^ and 
obtain the time-reversal symmetric open quantum map 

[HEl 



(3) 



which describes how the internal wave function tp is de- 
graded sequentially in successive decay attempts (where 
the integer n is the stroboscopic time). 

If we consider a short excitation pulse which leaves 
the system in the internal initial state \a) = V'(0)i the 
quantum map ^ yields 
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for the forward and reversed time evolution of the system. 
By contrast, for a continuous excitation process, Eq. Q 
implies 



a{e) = 1 + 2Re (a| FinQF, 
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Both types of experiments are related to one-another by 
the fact that the Fourier transform of a{e) is just the 
return amplitude. That is, in the energy and time units 
of Eq. (dJ: a{t) = M(a|V'(n)), where t = n/M. The 
stroboscopic form factor is then obtained from Eq.(^. 
Note that the spectrum of F is homogenous on the unit 
circle, so that we may set L = M. 

Direct processes. We now incorporate quasi- 
deterministic direct processes by assuming that part of 
the initial amplitude escapes during the first iteration of 
the map. Note that the total cross section contains the 
forward, but also the backward (time reversed) evolution, 
Eq. This means that the autocorrelation function 

will be equally sensitive to direct decay in either direc- 
tion. To study this dependence quantitatively, we assume 
the following decomposition for the initial state: 



|a) - \ao) +F^^^p\a+)+Fi,,p\a-) 
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Here, |ao) and \a^) are assumed to be real in the basis 
chosen. Also, we assume that the three terms are approx- 
imately orthogonal to each other, e.g. {a'^\p^ Fp\a~) « 
which would follow from the absence of any prompt pro- 
cesses in the full scattering svstem. |30j Hence, we find: 
||a||2 = ||ao||2 + ||c,-||2 + ||„+||2^ While |ao) leads to 

purely indirect decay, ||q;'''|P gives the probability for di- 
rect decay within the first step of the open map, Eq. JSJ, 
forwards/backwards in time. If ja"*") = |a~), \a) becomes 
real, and the whole decay process becomes symmetric in 



time. The orthogonality of the three terms in Eq. © 
then implies that the maximum amount of direct decay 
(in forward or backward evolution) is restricted to one 
half. 

To make contact with RMT we assume that Fout is 
taken from the circular unitary ensemble, which then 
implies that F = Fout-F^^ is a member of the circu- 
lar orthogonal ensemble The resulting scattering 
ensemble [i^l is essentially equivalent to the random 
Hamiltonian ensemble discussed in 0. In particular, the 
S-matrix correlations are well described by the VWZ- 
integral |2^ . This allows to calculate the autocorre- 
lation function of the cross section in the mapping for- 
malism of Refs. The result is (for details see 

H) 

C{t) = \\ao\\'Co{t) 

+ l\\aor{\\a+r + \\a-r + 2{a+\a~)) Coiit) 

+ \{\\a+r\\a-r + {a+\a-r) C,^it) , (7) 

where Coi(t) and Cii(i) are given by the expression jSJ, 
but with the function /q replaced by /oi ~ t and /n = 
{x + r)/{l + 2r + x) + {t-r)/{l + r-t), respectively. The 
auxiliary correlation functions Co{t), Coi(t) and Cii(t) 
all tend to two as f — > 0. The functions Co{t) and Cii(t) 
can be seen in Fig. ^ They correspond to the limit cases 
of purely indirect (||q;*|P = 0) and maximally direct 
(||a±f = 0.5) decay. The function Coi(t), if plotted, 
would lie in between. On a time scale given by tn 3> tEhr, 
all three functions are notably different. Equation ^ is 
the central result of this paper. It generalizes related 
results obtained within the autonomous random-matrix 
model 0, 01 . These results are recovered (for ballis- 
tic decay) by choosing ja"*") = |q;~), while the universal 
random-matrix prediction J^jl is obtained for |a*) = 0. 
In general, the form factor depends on three indepen- 
dent parameters, the probabilities for direct decay 
in forward/backward time evolution as well as the angle 
between the two corresponding amplitude vectors. 

Maximally asymmetric decay is obtained by setting 
\a~) to zero. It yields: 

Cit) = WaoW Coit) + \ ||aoini«+f C^i{t) . (8) 

This expression also describes the decay form factor for 
an initial state \a) ~ Fin (|ao) +p\oi^))i and yields the 
same dynamics as choosing the initial state |ao) +P |q^~) 
for the open map x(n -1-1) = Q F x('T') used e.g. in 

Rcf. [mil. 

Numerical results. The decay form factor of realistic 
physical systems can be obtained from experiments or in 
sophisticated numerical computations [22|. Here we test 
our theoretical predictions by numerical investigations of 
a simple model system, the open kicked rotator [l^ll3| . 
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FIG. 1: (color online) Form factor {C{t)) times exp{Nt) for 
the case of time-reversal invariant decay (real initial coher- 
ent wave-packet in the open kicked rotor) and for different 
amounts of overlap of the initial state with the regions of 
fast decay (solid lines). Here, M = 4096 and iV = 50. The 
corresponding theoretical results [Eq.JTJ, with \a'^) — \a~)] 
are shown with dashed lines. From top to bottom ||a:*||'^ ~ 
0,0.25,0.417,0.5. 

which was used in most studies of the stationary prop- 
erties of quantum systems with a finite Elircnfcst time 
[M ^3 ■ In the position basis the Floquct matrix can 
be decomposed as F — X W H U X , where 

n = diag(e-*'^^-'/^^), 

X = diag(exp[-iA/V"(27rfc/A/)]), (9) 

and V{e) = K{cose - 7 sin 20) with 7 = 1. For 7 = 0, 
the standard kicked rotator, there are three different sym- 
metries present: time reversal, reflection and conjugation 
(the combination of time reversal and reflection) . A finite 
value of 7 breaks the reflection and also the conjugation 
symmetry. The initial wave-packet |a) is constructed as 
a superposition of two coherent states, located around 
the line p = and therefore real. The amount of overlap 
with the regions of fast decay (equal for both backward 
and forward evolution) was controlled by changing the 
location of and/or squeezing the coherent states. 

We analyze the normalized autocorrelation function 
C{t) = 2 (7[cr](t)/C[cr](0), whose hmit value at smaU 
times is two, independent of the probabilities for direct 
decay. This makes sure that the asserted sensitivity to 
direct processes is related to the shape of the autocorrela- 
tion function, not its norm (in practice, the cross sections 
measured or calculated usually lack a reliable absolute 
scale). 

We first consider the case of a real-valued initial state, 
with symmetric decay in time. Figure ^ shows the re- 
sults for the form factor {C{t)), scaled by the classical 
survival probability exp{—Nt). The numerical results are 
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FIG. 2: Numerical data and theory for Y{t), for N = 20 
channels (left panels), and A'' = 50 channels (right panels). 
The overlap of the initial state with the regions of fast decay is 
given by \\a~f = 0(a), 0.595(b), 0.907(c), 0(d), 0.571(e), and 
0.906(f). Areas in gray represent the range of the standard 
deviation of Y{t), while areas in black represent the standard 
deviation for {Y{t)) when the average is taken over 100 sam- 
ples. The dashed lines show the theoretical curves for the 
different values for ||q~|P; the respective theoretical curve 
corresponding to the numerical result is plotted in light gray. 



obtained for a large number of internal modes M ~ 4096 
(required for our random-matrix model to apply) and a 
large number of open channels = 50. This allowed us 
to reach = 0.417 as the maximal probability for 

direct decay. We take averages over 100 different real- 
izations of the system, each corresponding to a different 
value of the kick-strength K (for all of them the classical 
dynamics is completely chaotic). The theoretical results 
derived from Eq. |(7J) with parameters corresponding to 
the numerical data are shown with dashed lines. The nu- 
merical results agree perfectly well with the theoretical 
predictions, within the remaining statistical uncertainty. 

So far we have only discussed an ensemble average, 
while experimentally or numerically it is not always pos- 
sible or desirable to change any system parameters. Can 
one detect the direct quasi-deterministic decay for an in- 
dividual system? This question is addressed, for the case 
of maximally asymmetric decay, in Fig. [3 We show the 
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results for M = 10000 and two values of A'': N = 20 in 
panels (a), (b), (c) and = 50 in panels (d), (e) and (f). 
On the vertical axis we plot the function Y{t) which we 
obtain from the numerical data through the transforma- 
tion m 



Y{t) 



Nt 



C{m/M) 



* t^o(C^o(™/M)-C^oi(WM)) 



which allows to obtain a best fit value for 



(10) 



by linear 



regression (not shown). In this sum we have discarded 
the first 20 kicks to exclude system-specific processes at 
the Ehrcnfcst time scale. It removes the fluctuations in 
time, so that only the samplc-to-sample fluctuations sur- 
vive. These arc shown by the gray areas in Fig. 2, while 
the expected standard deviation for the average over 100 
samples is shown in black. The dashed lines show the dif- 
ferent theoretical curves, obtained from Eq. (jSJ, that fit 
best each one of the numerical results. Since the sample- 
to-sample fluctuations (light gray) often do not cover the 
whole available range between the theory for purely in- 
direct and dominantly direct decay, even a single exper- 
iment may be sufficient to extract some information on 
the presence of direct decay processes. 

We have developed a dynamical description for the 
form factor C{t) of a decaying quantum system, taking 
into account the effect of direct (sub-Ehrenfcst-time) de- 
cay processes. We derived a general analytical expression 
for C{t) at times of the order of the Hcisenbcrg time, on 
the basis of a suitable random-matrix model. While ear- 
lier studies 0, 0, 0, 0, 0| implicitly assumed that 
the wave packet will follow the same dynamics in for- 
ward and backward evolution, we allow for excitation 
mechanisms which lead to asymmetric decay in time. In 
physical terms, asymmetric decay would result from the 
preparation of an initial wave packet with finite group ve- 
locity. Properly directed this might strongly enhance the 
probability for prompt decay, but only in one direction 
in time. 
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